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Abstract 

We consider Anosov thermostats on a closed surface and the X-ray trans- 
form on functions which are up to degree two in the velocities. We show that 
the subspace where the X-ray transform fails to be s-injective is finite dimen- 
sional. Furthermore, if the surface is negatively curved and the thermostat is 
pure Gaussian (i.e. no magnetic field is present), then the X-ray transform is 
s-injective. 

1 Introduction 

Let M be a surface, a smooth, closed, orientable 2-manifold with a Riemannian 
metric, and let </>t : SM — > SM be a flow on its unit sphere bundle. 

Given a smooth function g : SM — > R, the X-ray transform of g is the map 
determined by 1(g) (T) = j T g, where V is any closed trajectory of (f>. Is it possible 
to recover g if we just know its integrals along every closed orbit? Since / is linear, 
this is obviously equivalent to asking when / is injective. In general / is not injective, 
because non-zero functions of the form g = d/dt(u o <f>t)\t—0 (called coboundaries) , 
where u : SM — > K is a smooth function, lie in the kernel of the X-ray transform. 
In fact, if the flow tp t is Anosov, the smooth Livsic Theorem [8] tells us that g is in 
the kernel of the X-ray transform if and only if g is a coboundary. 

Classically, one considered recovering the function g in the special case g S 
C°°(M), a function of position only rather than both position and velocity. In the 
case g is linear in velocity we can write g(x,v) — gi(x)v'\ and the study of the 
X-ray transform has applications to Doppler tomography, photoelasticity and other 
physical problems [15] . 

In this paper we are interested in studying injectivity properties of the X-ray 
transform for functions g which are the restriction to SM of a symmetric 2-tensor 
q plus a 1-form a. The motivation for looking at the kernel of / for this class of 
functions comes from rigidity questions in dynamical systems and spectral geometry. 
Symmetric 2-tensors and 1-forms appear naturally when one linearises rigidity ques- 
tions involving Riemannian metrics and magnetic fields (cf. and [5]). There is a 
similar circle of questions for manifolds with boundary (motivated by the boundary 
rigidity problem), but we do not consider these here. 

Quite a lot is known about injectivity properties of the X-ray transform for 
geodesic flows and in this paper we wish to extend these to more general flows on 
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unit tangent bundles. We consider thermostats; so for i a rotation by ir/2 in the 
tangent planes, a generic path 7 in M is determined by 

= Kl: 7) Hj where A(x, t>) = /(at) + (e(a:), w), (1) 

/ some function and e some vector field. To simplify the exposition we shall assume 
that e has zero divergence (this is no restriction at all, cf. [T2J, Lemma 2.2]). 

As above, <f> will denote the induced flow, and a simple calculation confirms the 
sphere bundle is invariant. Let F be the vector field induced by the flow. Notice 
that if / = and e = then the system reduces to the geodesic flow; if only e = 
then the system is a magnetic flow. 

Given a pair [g, cr], where q is a symmetric 2-tensor and a is a 1-form, we will 
say that [q, a] is a potential pair if there exist a function h on M and 1-form tp such 
that F(h + = q + a (we consider restrictions of all the objects to SM). We note 
that because of the special form of A in fTJ, F(h + ip) is of degree at most two in 
the velocities. Clearly any potential pair gives rise to a function in the kernel of /. 
In Section 4 we will show that for pairs [q, a] there is an orthogonal decomposition 
into potential and solenoidal pairs. We say that I is s-injective if I[q, a] = implies 
that [q, a] is a potential pair. Equivalently, the restriction of / to solenoidal pairs 
is injective. 

Here is the main question we wish to address: 
Question. If the thermostat <fi is Anosov, is I always s-injective? 

Note that if I[q, a) = 0, the smooth Livsic Theorem [5] gives a smooth function 
u : SM — > R such that F(u) = q + a and the main point we are discussing in this 
paper is whether u is at most of degree one in the velocities. 

For geodesic flows, V. Guillcmin and D. Kazhdan [10] proved s-injectivity when 
M has negative curvature and the extension to higher dimensions was done by C. 
Croke and V. Sharafutdinov [2]. A key ingredient in the latter proof was the Pestov 
Identity. Relaxing the constraint on the curvature, N. Dairbekov and Sharafutdinov 
showed that by using the Anosov structure one could modify the Pestov identity 
and still obtain positive results [7]. Here, s-injectivity was proved up to a finite 
dimensional subspace. 

For thermostats, Dairbekov and the second author [3] showed that the integrals 
along closed orbits determined a function r(x) + o~ x (v), where a is a 1-form, up 
to an exact 1-form; for functions of degree at most one in velocity, exact 1-forms 
form the space of potential functions. This result allowed the authors to show that 
an Anosov thermostat on a surface has zero entropy production if and only if the 
external field e has a global potential. 

Let P be the set of all smooth potential pairs and let Z be the kernel of / acting 
on smooth functions of the form [q, a] . 

Theorem A. For an Anosov thermostat on a surface, the inclusion P C Z has 
finite codimension. 

In other words I is s-injective up to a finite dimensional subspace. It is unknown 
if P = Z always. 
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A pure ( Gaussian) thermostat is one with no magnetic term. Thus 7 = A(7, 7) ij, 
where X(x,v) = (e(x),iv). Pure thermostats are reversible, like the geodesic flow, 
whereas magnetic flows are not. 

Theorem B. For a pure Gaussian thermostat on a negatively curved surface, the 
X-ray transform is s-injective. 

We note that thanks to a theorem M. Wojtkowski [19] a Gaussian thermostat on 
a surface of negative curvature is always Anosov (recall that we are assuming that 
e has zero divergence). However these thermostats could have very large geodesic 
curvatures (i.e. e could be arbitrarily large) and the usual direct approach using the 
Pestov identity does not work. Indeed, the bound needed is derived in [H Theorem 
4.6]; 

K-H{\) + \\ 2 <0, 
5 

where H is a vector field on SM defined in Section ® This is clearly not true 
for arbitrary A. We bypass this problem using the ideas of Sharafutdinov and G. 
Uhlmann in [17j . 

There is a key difficulty in extending Theorem A to higher dimensions. It is 
unknown if the weak bundles are transversal to the vertical fibration for thermostats 
in dimension > 3. This is essential in the proof of Theorem A and it is proved in 
dimension 2 in [3]. 

For magnetic flows this property is known |15j and in principle there should 
be no serious obstacles proving Theorem A for them using the Pestov identities 
developed in [5]. However, even the following remains open: 

Question. Consider an Anosov magnetic flow on a surface, so that e = in the 
equation of motion ([I]), and assume K — H(\) + A 2 < 0. Is the X-ray transform for 
this flow s-injective on pairs [q, a]? 

2 An integral identity 

Let M be an oriented Riemannian surface, so that the linear map i : T X M — > T X M 
is well defined as follows; for v £ T X M with unit norm we require {v,iv} to be 
an oriented orthonormal basis. Recall from basic surface geometry that the sphere 
bundle 7r : SM — > M is a principle circle bundle; see [18]. Let V be the infinitesimal 
generator of the circle action. 

Let X be the vector field over SM induced by the geodesic flow. Define a third 
vector field H using a flow <p orthogonal to the geodesic flow: given (x, v) S SM, let 
7(x,v ± ) (t) De the unique geodesic from x heading in a direction (with the orientation) 
perpendicular from v. Define <pt(x,v) — ( i y^ XtV ±\{t),u{t)), where u is the parallel 
translate of v along 7, and H the vector field on SM that ip induces. Notice that 

{d-K( x , v ){H),iv) = 1. 

A basic result in Riemannian Geometry says that T( x v jSM is spanned by 
{X, H, V}. Again referring to [18], we recall that 

[V,X]=H, [H,V]=X, [X,H]=KV, (2) 
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where K is the Gaussian curvature of M. 

As in [17], for real functions A and c on SM we define a modified basis {F, H c , V} 
of T^ X V )SM in order to explore flows other than the geodesic flow, with 

F := X + XV, H c := H + cV. (3) 

For example, if X(x, v) = f(x) then the flow generated by F is a magnetic flow, 
whereas if X(x, v) — lo x (v), for some 1-form to, then we have a Gaussian thermostat 
flow. The importance of c will become clear in Section [3] It is worth pointing out 
that for a path 7 induced by the vector field F, 

as in ([1]) in the introduction. 

We will need a curvature type function on the sphere bundle, K(x, v) := K(x) — 
H C (X) + X 2 . 

Remark 2.1. When X(x,v) — f(x), and so the flow is purely magnetic, this ex- 
pression reduces to the magnetic curvature of [T] if c = 0. 

Lemma 2.2 (Pestov identity). For a twice differ entiable function u : SM — > R, it 
is always true that 

2H c u ■ VFu = (Fu) 2 + (H c u) 2 - (F(c) + c 2 + K){Vu) 2 

where these terms will ( +F(H c u ■ Vu) + V(X)H c u ■ Vu ^\ 
be made to vanish later < —H c (Fu ■ Vu) — V(c)Fu ■ Vu 
after integration. [ +V(H c u ■ Fu), 

Remark 2.3. This lemma was given in [17] for the geodesic flow, A — 0. Without 
the modified horizontal derivative consideration, so that c = 0, a version of this 
lemma was given in [3]. 

Proof. Combining equations ([2]) and ([3]) we see that, with respect to the modified 
basis, the commutation relations are 

[V,F]=H c + (V(X)-c)V, 

[V,H c ] = -F+(V(c) + X)V, 

[F, H c ] = -XF - cH c + (F(c) + c 2 + K)V. 

A proof can be derived from an application of these relations. There is no difficult 
step once we know of the existence of such an identity, so the explicit algebra has 
been relegated to the Appendix. □ 

On SM we consider the standard Liouville volume form with associated Liouville 
measure fi. Integrating the Pestov identity (fj]) over SM with respect to fi, yields 
the key equality. The last five terms disappear because each flow induced by the 
unperturbed vector fields, X, H and V, leaves the volume form invariant. Thus for 
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any function / on SM we have J X(f)d[i = J fCxd^L = 0, and similarly for H and 
V. 

This argument immediately implies the integral of the last term in ^ will be 
zero. Integration by parts after applying the above argument can be used on the 
other terms. Using F = X + XV, for instance, 

/ F(H c u-Vu)d^L = {)+ I \V(H c u-Vu)dfi= - I V{\)H c u ■ Vud/j,. 

JSM JSM JSM 

Explicitly then, the key integral identity derived for use in Section [4] is 
2/ H c u-VFudfi=\\Fu\\ 2 L2 + \\H c u\\ 2 L 2 - [ (F(c) + c 2 + K){Vuf dp. (5) 

JSM JSM 

Now we can follow the exposition in [7J, where a function c : SM — > K was 
chosen so as to make the contribution of the last term in the above expression 
arbitrarily small. This was equivalent to using an approximate solution of the 
generalised Riccati equation. When the base manifold is a surface one can use an 
exact solution. In the next section we outline how and why this can be achieved. 



3 Anosov splittings and the Riccati equation 

In this section </> is the thermostat flow in the sphere bundle. We consider variations 
of a path 7(i) = iro<f> t (x, v), for some (x, v) G SM. Notice 7 is defined on M, rather 
than on SM as elsewhere in the paper. We let f(s, t) = ir(4> t (Z(s))) be a variation of 
7; / depends on a choice of curve in TM called Z, with Z(Q) — £ G T/ XjV \SM. The 
image of / is a small rectangle covering 7. We define the Jacobi fields (dependent 
on as J c = df/ds(0, ■). Then 

J € + i?( 7 , J 5 ) 7 - d\(J 5 , j 6 ) ij - = 0. (6) 

The analysis is relegated to Subsection IA.2I in the Appendix. 

The importance of Jacobi fields for this paper is contained in the next Lemma. 
We will need a few basic considerations of the tangent bundle, see [13] . 

On any manifold there is always a canonical vertical bundle V(x,v) := ker(c?7r( :E v \ ) 
over the tangent bundle. The choice of a horizontal bundle TL{x, v) - in order that 
H © V = TTM - corresponds to a choice of connection on TM. 

We use the Levi-Civita connection to define K. : TTM — > TM, the connection 
map. This will induce a choice of TC. Given £ e T( x v yTM , choose a path Z(t) = 
(a(t), z{t)) C TM adapted to f ; this means a(0) = x, z(0) = v and Z(0) = £. Then 

/C(0 := (V d *)(Q). (7) 
We set H(x, v) := kerfX^ „)). It is an easy exercise to check that 

T (x , v) TM -> T X M © T X M 
£ 1 ► (dirt 

is an isomorphism; throughout this section we use this splitting. It is also immediate 
that for the unperturbed vectors from the last section V G V while X,H G H. 
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Lemma 3.1. If we choose the curve Z C TM adapted to £, such that Z(0) = £ } 
and let be the Jacobi field varying with Z , then 



#*(£) = (Jf(t),J e (i)). 

Proof. 

dn o d4 t {£) = d(ir o t )(£) = ^-(tt o fa o Z(s)) = J e (t) 

d 



=o dt ds 



ds dt 



□ 



By definition, for an Anosov flow cf> t : SM — > there exists a strong stable 
bundle E~ C T(SM) and a corresponding constant < fi s < 1 such that 

WdUOW < lieile^S V£e£ s , Vi >0. 

Similarly, there is a strong unstable bundle E + , and T(SM) = E~ ® E a ® E + is 
a continuous splitting, where -E is generated by the flow direction. By a result of 
Hirsch, Shub and Pugh, p], the two weak bundles, E° © E~ and E° ® E+ , are of 
class C 1 . It is these weak bundles that we use in the construction. 

Pick a weak stable or weak unstable bundle and denote it E. We know that E 
is invariant under the flow. By the work of Dairbekov and the second author in [5], 
E is always transverse to V for a thermostat flow. Hence there exists a pointwise 
defined linear map S( XtV ) : TL — > V such that 

E {X , V ) = gr&ph(S {X:V) ) = (u, S( x>v) (u)) cHffiV~ T X M © T X M. 

Using the identification above, we view Si XtV ) as a map T X M — > T X M. 

Notice that as -E is a weak bundle, certainly it contains the vector field induced 
by the flow; E 3 F = X + XV. Thus S(dir(X)) = K(XV), i.e. S {x . v) (v) = Xiv. We 
define a new function on the sphere bundle which together with A determines the 
weak bundle, 

c(x,v) := (S( x>v )(iv),iv). 

Remark 3.2. Of course, this actually defines two functions; one for the weak stable 
bundle, and one for the weak unstable bundle. There should be no confusion as the 
arguments follow identically, and the brevity helps give a clear account. 

If we take £ in a weak bundle E, then d<f>t (£) S E Vt, as E is flow invariant. Since 
E can be represented as the graph of the bundle map 5 and d<p t (0 = (J^(t), J$(t)) C 
E by Lemma |3. 11 we have 
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We now proceed as usual [T3] and from ^ obtain 



F(c) + A 2 + c 2 + K - H c (\) = F(c) + c 2 + K = 0. 



(8) 



The explicit calculation is done in the Appendix. 

Remark 3.3. Since the weak bundles are C 1 , we do not need to approximate 
c in the derivation of the Pestov identity, as was required in [7] . This means that 
the Pestov identity on surfaces greatly simplifies when modified to any Anosov flow; 
see the next Lemma. In general, this will not hold in higher dimensions. 

Lemma 3.4. Let u : SM — > R be any function. From the discussion in this Section 
we can pick c : SM — > K ; corresponding to either the weak unstable bundle or the 
weak stable bundle, such that 



4 The energy estimates method 

Throughout, u and p are functions on the sphere bundle. In [2], Croke and Shara- 
futdinov studied the geodesic cohomological equation, 



on a compact negatively curved manifold. For p a symmetric tensor function of 
the form p(x,v) = Pi 1 ...i m .(x)v tl ■ ■ -v %m , it was shown that p must be potential. In 
that paper, a potential function was of the form p = dv, where d = d s = a o V 
was covariant differentiation of a tensor followed by symmetrisation of the indices. 
We now outline why we cannot expect such a simple constraint for solutions of the 
thermostat cohomological equation 



with F = X + XV, as in ([T]). We can uniquely split a function p on SM into its 
even and odd parts; 



such that p — p e +p°. For the general thermostats in this section, A has both even 
and odd components. By considering the even and odd parts of the thermostat 
cohomological equation we immediately realise the extra XV term requires a more 
involved definition of the operator d. Similarly, we should not restrict our attention 
just to symmetric tensor fields of some fixed degree. 

We are interested in the X-ray transform / applied to functions p of the form 
p{x, v) — q x (v, v) + a x (v ), where q is a symmetric 2-tensor, a is a 1-form on M . If D 
is a notation for a function space (C k , LP, H k , etc.), then we will denote by D(M) 




Xu = p, 



Fu = p, 



2p e (x, v) =p(x, v) + p(x, —v) 
2p°(x,v) =p(x,v) ~p(x, -v), 
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the corresponding space of pairs f = [q, a] , with q a symmetric covariant 2-tcnsor 
and g a 1-form, and denote by V{M) the corresponding space of pairs w = [ip, h], 
with -0 a 1-form and h a function on M. In particular, L 2 (M) is the space of square 
integrable pairs f = [q, a], and we endow this space with the norm 

||f|| 2 = / {\q\ 2 + \a\ 2 }dvol, (9) 

JM 

with the corresponding inner product. In the space C 2 (M) we will consider the 
norm 

||w|| 2 = I (\^ 2 + h 2 ) dvol. (10) 

JM 

Clearly, the norm of a pair [q,<r] in L 2 (M) is equivalent to the norm of the 
corresponding function p(x,v) in L 2 (SM). We wish to construct an operator d 
such that I(d[tp,h]) = for all [ip,h] £ H 1 (M). For this it is enough to have 
d[ip, h] = [q, a], where q + a = F(h + -0). 

We use the hypothesis X(x, v) = f(x) + (e, iv); earlier in the paper A could have 
been any once differentiable function on SM. Set 8 x (v) := (e(x),v), so A = f+V(8). 

We compute (X + \V){ip + h) = d s i; + dh + f ■ V(ip) + V (9) ■ V(ip) + 0. Projecting 
onto the odd and even parts, this is rewritten as 

A/A _ fd s i; + V(9) ■ Vtp\ _ fd s + V{6) -V \ (ip\ 

d {h) -{ dh + f-V^P ) ~\ f-V d-v){h)- 

Hence d is the operator 

(d* + V{6) -V \ 

\ f-v d-vj- 

Definition 4.1. We call a pair [q,a] £ L 2 (M) potential if the equations d[tp,h] — 
[q, a] hold with [ip, h] £ H 1 (M). A pair [q, a] is called solenoidalif it is L 2 -orthogonal 
to all potential pairs. 

We let { • } s denote symmetrisation. If a pair [q, a] is orthogonal to all potential 
pairs, then for all tp an d h we have: 

/ (q, d s ^ + {V(6) ■ V(iP))} s + {a, f ■ V{^) + dh) dvol = 0. 

JM 

Then 

- f (6q + V{i le q) + / • V(a), + {5a, h) dvol = 0, 

JM 

where 8 is divergence and ii e q is the 1-form given by (ii e q)x( v ) — qx(^ e ( x ),v). 
Therefore, 



6q + V(t ie q) + f-V(a)=0, 
(5cr = 0. 
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Thus [g, a] is solenoidal if S[q, a] — 0, where 

(5 + Voi ie f-V 
^ 5 

and d = — S* . The operator — 8d acting on pairs w = [ip, h] is elliptic (the proof is 
the same as in [5J Section 3.3] for the pure magnetic case) and if the thermostat <p 
is transitive, its kernel is given by the pairs [0, h] where h is a constant. Indeed the 
kernel is given by those smooth w such that dw = 0, which means that h + ip is a 
first integral of the thermostat. Assuming transitivity we must have ip = and h 
constant. 

Suppose now that f = [q, a] is given. By solving the equation ddw = Sf for w 
and setting f s := f dw we see that we can decompose f as 

f = dw + f s , Sf s = (11) 

where f s is uniquely determined and w is uniquely determined up to pairs of the 
form [0, h) with h constant. If f is smooth, then f s and w are also smooth. Finally 
we note that an Anosov thermostat on a surface is always transitive by a result of 
E. Ghys 0. 

To simplify notation, we let |u| denote the L 2 norm of u, and ||w|| its H 1 norm. 
Proposition 4.2. If u G C°°(SM) satisfies the cohomological equation 

p = Fu, 

with pix, v) = q x (v, v) + a x (v), then there exists a constant C such that 

\\uf <c(\u\ 2 + \S[q,a}\ 



Proof. Substituting into the Pestov integral identity of Lcmma l3.4l (with c a solution 
of the Riccati equation as found in Section [3]) , we see 

H c u -Vpdfi= \p\ 2 + \H c u\ 2 . 

SM 

The initial steps are algebraic manipulation. All integrals are over SM with respect 
to the Liouville measure /z. 

\p\ 2 + \H c u\ 2 = 2 J H c u-Vp 

= 2 J H c {u ■ Vp) - 2 J u ■ H c V(p) 

= 2 J cV(u ■Vp)-2Ju- {HVp + cV 2 p) 

= -2 J u-(Vc-Vp + HVp + cV 2 p), 
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\p\ 2 + \H c u\ 2 < - 



2 J u- (HVp + Vc-Vp + cV 2 p). (12) 



We look to bound each term on the right hand side. Note that since SM is compact, 
the continuous function c : SM — * R is bounded. Since it is of class C , Vc is also 
a bounded function. 

We have a prescribed form for p, which gives 

(Vp)(x,v) = V(q x (v,v) + a x (v)) = 2q x (iv,v) + a x (iv). (13) 

So \Vp\ < 2\p\, and similarly \V 2 p\ < 4|p|. 

Let us bound the integral of uHVp. Using (fT3")) there are two integrals to 
consider, that of uHVa and that of uHVq. We will need the following easy results 
(cf. [HI Lemmata 5.2 and 6.1]); 

X(a) + HV(a) = 8a, X(q) + HV(q)/2 = Sq. (14) 



J u ■ HV{a) = J u-Scr- J u- X(a) 
u ■ 5<t + / X(u) ■ a 



u ■ 5a + j p ■ a — j XVu ■ a 
u-(5a + V(X)a + XV(a)) + j p ■ a 
u ■ (6a + 2/ • V(a)) + [ u ■ V(V{6) ■ a - f ■ a) + / p ■ a. 



J u ■ HV(q)/2 = J u- Sq~ J u- X{q) 

= j u ■ Sq + J p ■ q — J XVu ■ q 

u-(Sq + V{i ie q)) + / u ■ (fV(q) - L e q) + I p ■ q. 



Notice that the first integral in each of the above calculations - switching back 
to the original coefficients in (fT2|) - will sum to yield 



-2 Ju-(Sa + 2f-V(a))-A J u ■ (Sq + V{i ie q)) < A\u\ ■ \5[q, a]\. (15) 
The last term in each of the two expressions are also bounded, 
-2 J p-a-A J p.q = -4 J a 2 /2 + q 2 < 0. 
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Finally, the middle integral of both expressions is bounded by C'\u\ ■ \p\, by the 
same arguments as when dealing with the operator V in (|13p . 
Feeding this into (fT2")) we have a bound on the left hand side: 

\p\ 2 + \H c u\ 2 < -2 J u-{Vc-Vp- cV 2 p) - 2 J u ■ HVp 

< C\u\ ■ \p\ + A\u\ ■ \S[q, a}\+0 + C'\u\ ■ \p\. 

Here C and C are constants, dependent only on A and c. Hence the constants are 
determined by the choice of Anosov flow, rather than u and p. Let C denote some 
generic constant greater than 1. 

\p\ 2 + \H c u\ 2 <c(\u\-\p\ + \u\-\8[q,a}\). 

Let c u and c s be the two Riccati solutions from Section [3J corresponding to the 
weak unstable and stable bundles of the Anosov flow. Adding their corresponding 
inequalities when substituted into the last inequality, we have 

M 2 + \p\ 2 + \H c uu\ 2 + \H c su\ 2 < C(\u\ 2 + \u\ ■ \p\ + \u\ ■ |%,ct]|) . (16) 

We claim that the first Sobolev norm of u is smaller than the left hand side. Recall 
that ||u|| 2 := \\u\\ 2 H1 = \u\ 2 + \Xu\ 2 + \Hu\ 2 + \Vu\ 2 . Using the definition of F and 
the cohomological equation, 

|A w | 2 <b| 2 + |A^| 2 <C(H 2 + H 2 ). 

Since the flow is Anosov, c s ^ c u everywhere, which means the linear system 

(Si) ■ C (?) 

has a unique pointwise solution for H and V . Say H = aH c u + [3H c s . Then 

\Hu\ 2 < a 2 \H c uu\ 2 + f3 2 \H cS u\ 2 < C (\H c u U \ 2 + \H c s U \ 2 ) . 

The bound for \Vu\ 2 follows in the same way. 
Therefore, using (TT6"]) . 

\\u\\ 2 <c(\u\ 2 + \u\-\p\ + \u\.\S[ q ,a}\). (17) 

Finally, we use the cohomological equation to bound \p\ = \Fu\ < C\\u\\ which 
turns (|17p into a quadratic inequality in ||u||. Solving gives the required bound. □ 

We now have the same bound on the first Sobolev space as in [7] , but for Anosov 
thermostat flows on a surface. Thus we have the same consequence as their paper. 
We quickly outline the proof they gave. 

Theorem A. For an Anosov thermostat on a surface, the inclusion P C Z has 
finite codimension. 
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Proof. If the inclusion is not finite dimensional, take a linearly independent sequence 
of functions {zk}^ = i 6 Z\P. Using the potential-solenoidal decomposition (fTTjl . for 
each k there exists functions such that 

z k = dx k +y k , Sy k = 0. 

But z k £ Z yk S as dxfc is potential and therefore integrates to zero along 
all trajectories. Using the Livsic Theorem [5], we can find another sequence of 
functions {u k } k x L 1 with 

Fu k = y k . 

Using linearity and Gram-Schmidt, we can assume {u k } k x L 1 is an orthonormal set. 
The y k may also need to be adjusted, but by linearity they will still be solenoidal. 

Applying the key inequality, Proposition 14. 21 on Fu k = y k the Sy k term disap- 
pears: 

IKH < chi = c, 

as the u k were normalised. Therefore in H 1 the set {u k } is bounded. This leads 
to the required contradiction, as the embedding H 1 <—* L 2 is compact, which would 
imply the infinite orthonormal sequence u k had a convergent subsequence. □ 

5 Pure thermostats in negative curvature 

For geodesic flows on surfaces, Sharafutdinov and Uhlmann J7] found a useful new 
form of the integrated Pestov identity ([5]) . The main advantage of their modification 
was that they had effectively 'completed the square', controlling the (Hu, VXu)ip. 
term. Thus the full power of the identity could once again be applied, as was done 
i" El- 

Here we exhibit a subset of thermostat flows for which the results of [T7] still 
hold, so that the inclusion P C Z is an equality, P — Z. We proceed on a negatively 
curved surface with a pure thermostat flow, such that A = V[ff). The equation of 
motion takes the form 

— = (7, ie) ij 
and we will assume that e is divergence free. 

Remark 5.1. Notice that the odd and even parts of the cohomological equation 
are uncoupled, unlike in Section [JJ 

F(u e )=p°, F(u°)=p e . 

Again, we are considering coboundaries which are up to quadratic in the velocities. 
However, thermostatic coboundaries of up to linear order were understood in [21 
Theorem B]. This leads us to consider coboundaries which are purely even, p = q, 
for some symmetric 2-tensor q. 
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Lemma 5.2. Suppose F(u) = p, where u 6 C°° (S M) , p is a symmetric 2-tensor 
and F = X + V(9) ■ V . Let ip = V 2 (u) + u. Then the integral Pestov identity {3P 
can be rewritten 

\F<p\ 2 + \H c p\ 2 + \2-H c p- (9 + c)-Vp\ 2 = f H c -(V<p) 2 d» (18) 

JSM 

where, as before, | • | denotes the L 2 norm and we have a Riccati like function 
K c = F{6 + c) + (9 + 2c) (6 + c) + K. 

Proof. 

VFu = FVu + [V, F]u= FVu + H c u -(6 + c)Vu, 

and similarly 

V 2 Fu = FV 2 u + [V, F]Vu + H c Vu + [V, H c }u 
-(V(9) + V(c))Vu -(9 + c)V 2 u 
= FV 2 u + H c Vu -{9 + c)V 2 u + H c Vu -Fu+ (V(c) + V{9))Vu 
-(V(9) + V(c))Vu -(9 + c)V 2 u 
V 2 p = FV 2 u + 2H c Vu - 2(9 + c)V 2 u - p, 

where we have used Fu = p for the last line. 

Denote ip = V 2 u + u. Using the equation for V 2 p, 

Fip = FV 2 u + Fu 

= V 2 p - 2H c Vu + 2{9 + c)V 2 u +p + p 

VFp = -2VH c Vu + 2(V(9) + V(c))V 2 u + 2(9 + c)V 3 u + V 3 p + 2Vp 
= -2H c V 2 u + 2FVu + 2(9 + c)V 3 u + V 3 p + 2Vp 
= -2H c V 2 u + 2VFu - 2H c u + 2(9 + c)Vu 

+2(9 + c)V 3 u + V 3 p + 2Vp 
= -2H c ip + 2(9 + c)Vip + V 3 p + Wp. 

As mentioned in Remark 15.11 we have assumed a prescribed form for p ; that it 
is even and up to quadratic in the velocities. But then V 3 p + Wp = 0, and so 

VFp = -2H c tp + 2(9 + c)Vip. (19) 

But then 

2\ H c tp-VFipdn = -4:\\H c (p\\ 2 L2 +4 [ H c p- (9 + c)Vtpdn 

JSM '~ JSM 

= -pBjp -(9 + c)Vp\\ 2 L2 + [ (9 + c) 2 (Vp) 2 dfi 

JSM 
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Substituting this into the Pestov identity, (O in Section [2j transforms the left 
hand side when applied to the function <p = V 2 u + u: 

-\\2H c tp- (e + c)V(p\\ 2 L2 + [ {e + c) 2 {Vip) 2 d^ 

J SM 

= \\F V \\ 2 L2 + \\HM\h - I + c 2 + K){Vip? dpi, 

JSM 

and all that is left is to tidy the two integrands into one Riccati type function. 

F(c) + c 2 + K+(9 + c) 2 = F(c) + c 2 + K- H C V{6) + (V9) 2 + (9 + c) 2 

= F{9 + c) + {9 + 2c) (6> + c) + K, 

as e is divergence free, so we can use X9 + HV9 = div e = 0. □ 

Lemma 5.3. The two functions induced by the stable and unstable weak bundles of 
the flow, c s and c u of Section [5J are such that everywhere 

(9 + c s ){9 + c u ) < 

on a negatively curved surface. 

Proof. The fact that pure Gaussian thermostats are reversible gives c s (x, — v) — 
~c u {x,v). Thus 

9(x, -v) + c s (x, -v) = -(0(x, v) + c u {x, v)) 

and it suffices to show that the functions 9 + c s ' u do not vanish. We know that 

F{9 + c s ' u ) + c s ' u {c s > u + 9) + K = 0. 

Suppose 9 + c s vanishes at (x,v). Since K < 0, then F{9 + c s ) > at (x,v). It 
follows that {9 + c s )((f>t(x,v)) vanishes for at most one value of t G M. Recall that 
by a result of Ghys |9J the closed orbits of <j> are dense, hence we can find (y, w) 
close to (x, v) such that 9 + c s vanishes at (y, w) and the orbit of (y, w) is closed. If 
T is the period of the closed orbit we have (9 + c s )(y,w) = (9 + c s )(cj>T(y,w)) = 
which is a contradiction. 

□ 

Proposition 5.4. For a thermostat on a negatively curved surface, there exists a 
C 1 function c : SM — > K such that 1Z C is everywhere strictly less than zero. 

Proof. We already have two functions, c s and c u induced by the stable and unstable 
weak bundles of the flow, which satisfy a thermostat type Riccati equation ([8|). After 
a little rearrangement (and since e is solenoidal), and then adding the two almost 
identical equations we get 

F(c s +c u + 29) + c s (c s +9)+ c u {c u + 9) + 2K = 0. (20) 
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Recall K c = F(9 + c) + (8 + 2c) (0 + c) + K. Substituting 2c = c s + c u into flU}, 

= F(2c + 29) + c s (2c - c u + 9) + c u (2c - c s + 9) + 2K 
= F{2c + 26>) + (c s + c n )(2c + 9)- 2c s c u + 2K 
= F(2c + 29) + (2c + 20)(2c + 9)- 28(2c + 8)- 2c s c u + 2K 
= 2K c -2(8 + c s ){8 + c u ). 

Hence Ti-( c s +c u )/2 = (9 + c s )(8 + c u ), which is less than zero by the previous 
Lemma. □ 

Theorem B. For a pure Gaussian thermostat on a negatively curved surface, the 
X-ray transform is s-injective. 

Proof. Let Fu = p = q + a. Since the flow is Anosov, the simpler cohomological 
equation 

F[u) = a 

has a solution iff a = F(h) for some h £ C°°(M), using 3, Theorem B]. So we 
replace u by u — h, and assume p — q. 

Using Proposition 15.41 in (fT5)) of Lemma IB~2l with ip = V 2 u + u, 

\F^\ 2 + \HM 2 + |2 • H e <p -(6 + c)- V V \ 2 < 0. 

In particular, Ftp is zero and so ip is an integral of the flow. But the flow is Anosov 
which implies the only integrals are constants. 

We have V 2 u + u = k, some constant, which can be rewritten 

V 2 {u- k) + (u- k) = 0. 

For (x, v) £ SM we can write v — (cosw, sinw) 6 T X M for some real u>. Then there 
exists a, (3 6 C°° {M) who induce %j) = ip(a,P) £ Q 1 (M) such that 

u(x, v) — k + a(x) cos to + (3{x) sinu = k + ip(v). 

Thus u = k + ip, a constant plus a one form. Hence for the general coboundary 
p = q + cr, 

Fu = p => u{x, v) = h(x) + ipx{v). 

□ 

Remark 5.5. The hypotheses in Theorem B can be somewhat relaxed. All that is 
needed is that the thermostat cf> is Anosov and (9 + c s )(8 + c u ) < everywhere. 

Remark 5.6. There is an alternative approach to Theorem B. Let c := —9. For 
this choice of c, a simple calculation shows that 

[V, F] = H c , [V, H c ] = -F, [F, H D ] = -V(9)F + 9H C + (K - dive)V 
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and F(c) + c 2 + K = K — div e. If we assume e to be divergence free and we apply 
the Pestov identity ([5]) to if, we obtain 



2/ H e <p>VF<pdii=\\Fvf vt + \\H c <p\ ] (v- [ K(V<p) 2 d». 




But for our choice of c, VFip 



2H c tp (cf. HH)) and thus 



5||if cV ?||f 2 = ||f>]]f 2 - / K(V^ 2 dn 



JSM 



Hence if K < we see that F(cp) — H c (ip) = which implies that ip is constant. 
From this, Theorem B follows as above. 

A Appendix 

A.l A Pestov identity for thermostats on surfaces 

We modify the standard vector fields over SM, {X, H, V} as defined in [T5], to 



with A and c both functions on SM. Notice that F is the vector field induced by 
the thermostat flow 7 = A(7, 7) ij, and c could be chosen such that F and H c span 
a weak bundle if the bundle were always transverse to the vertical. 
By a trivial calculation, the modified commutation relations are 



where R = K - H c {\) + A 2 . 

Lemma A.l. Rewritten in terms of the basis {F, H C ,V}, the Pestov identity be- 



F = X + XV 
H c = H + cV 
V 



[V,F] 



H c + (V(A) - c)V, 



[V,H C ]-- 
[F,H C ] = -\F 



-F + (V(c)+X)V, 

cH c + (F(c) + c 2 + K)V, 



comes 



2H c u ■ VFu = {Fuf + (H c u) 2 - (F(c) + c 2 + K)(Vu) 2 
+F{H c u ■ Vu) + V(X)Fu ■ Vu 
-H c (Fu ■ Vu) - V(c)Fu ■ Vu + V(H c u ■ Fu). 
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Proof. 



2H c u ■ VFu - V(H c u ■ Fu) 
= H c u ■ VFu - VH c u ■ Fu 

= H c u ■ (FVu + [V, F]u) - Fu ■ (H c Vu + [V, H c ]u) 
= H c u ■ (FVu + H c u + (V(X) - c)Vu) 

-Fu ■ (H c Vu -Fu+(X + V{c)) ■ Vu) 
= (H c uf + (Fu) 2 + H c u ■ FVu + (V(\) - c) ■ Vu ■ H c u 

-Fu ■ H c Vu - (A + V(c)) -Vu-Fu 
= (H c u) 2 + (Fu) 2 + F(H c u ■ Vu) 

-FH c u ■ Vu + H c Fu ■ Vu - H c (Fu ■ Vu) 

+ ((V(X) - c)H c u - (X + V(c))Fu) ■ Vu 
= (H c u) 2 + (Fu) 2 + F(H c u ■ Vu) - H c (Fu ■ Vu) 

+([H C , F]u + (V(X) - c)H c u +(X + V(c))Fu) ■ Vu 
= (H c u) 2 + (Fu) 2 + F(H c u ■ Vu) - H c (Fu ■ Vu) 

- ((F(c) +c 2 + K)Vu - V(X)H c u - V(c)Fu) ■ Vu 
- (Fu) 2 + (H c u) 2 - (F(c) +c 2 + K)(Vu) 2 + 

F(H c u ■ Vu) + V(X)Fu ■ Vu 

-H c (Fu ■ Vu) - V(c)Fu ■ Vu + V(H c u ■ Fu). 



A. 2 Jacobi fields for thermostats 

Here <j) t is the thermostat flow in the tangent bundle. We consider variations of a 
path 7(t) = 7r o <fr t (x,v), for some (x, v) G SM. Notice 7 is defined on M, rather 
than on SM as elsewhere in the paper. 

We let f(s,t) = ir((f)t(Z(s))) be a variation of 7; / depends on a choice of curve 
in TM called Z, with Z(Q) = (x 7 v). The image of / is a small rectangle covering 
7. Further, we define nearby paths 7^ = f(s, ■) and the Jacobi fields (dependent on 



Using D/ds\ s =o — Vj on vectors over M and that A = X(j s (t),j s (t)), a function 



□ 



Z) as J = df/ds(0, ■). 

By the definition of curvature 



D Ddf 
ds dt dt 



DDdl d£ df df 
dt dtds [ dt' ds> dt 
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on SM, we have 



D_(Pis 



ds V dt 



) = —(\ii s )\ s= o = — A(7 a (t),7 a (t))| a =o -«7 + + AiVj(7) 
/ s=o as as 

Dd/(s,t) 



: dA( J, J) i7 + A i- 

dt ds 

dX(J,J) ij + XiJ. 



s=0 



On the first line we made use of the fact that the 1-1 tensor i is parallel with respect 
to the Riemmanian metric. Thus the generalised Jacobi equation is 

J + R(i, J)7 - d\(J, j) ij - XU = 0. 

Notice that dX(J, j) = (horizontal lift of J) (A) + (vertical lift of J) (A). 

As we are working on a surface we can study a general Jacobi field in terms of 
two real variables of t, say x and y, with 

J(t)=x(t)i(t) + y(t) » 7 (t). 

The horizontal lift of J is then xX + yH . Further, if the Jacobi field has initial 
condition £ G Ti xv \SM then, (J, 7) = which implies x = Xy and that the vertical 
lift of J is (y+Xx)V . Considering the generalised Jacobi equation in the ij direction 
then, 

y + Xx + X 2 y + Ky- xX(X) - yH(X) - (y + Xx)V(X) = 0, 

y~V(X)y + (K -H(X) + X 2 )y = 0, 

since X(X) + XV (X) — F(X) — A, and the analysis is complete on substituting 
c(t) = y/y. Note that in the notation of Section [3l 

S(J) = J = xS(j) + ySfyj) = (Xx + cy)ij = (Xx + y)«7. 
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